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Abstract 

We investigate a potentially large CP violating asymmetry in the decay of a 
neutral scalar or pseudoscalar Higgs boson into the tt pair. The source of the 
CP nonconservation is the complex mixing in the stop ti, r sector. One of the 
interesting consequence is the different rates of the Higgs boson decays into CP 
conjugate polarized states. 
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1 Introduction 



The standard model (SM) of particle interactions contains one CP violating parameter, 
which is a complex phase in the quark sector of the SM. This phase appearing in the 
quark mixing matrix of the charged current is expected to account for the observed CP 
violations in the K-K mixing, in the K decays, as well as in the potential CP violation 
in the B-B system. 

However, it is generally believed that new physics beyond the SM must exists. One 
of the major motivations for this is to understand the seemingly unnaturalness of the 
Higgs mass at the electroweak scale in the SM, the so-called gauge hierarchy problem. In 
addition, due to the difficulties of the SM to account for the baryon asymmetry of the 
universe as well as to resolve the strong CP problem, it is widely accepted that new sources 
of CP violation are needed. The most popular extension of the SM that addresses the 
hierarchy problem is the supersymmetric standard model|jI|, The extension has many 
more new (super-) particles and parameters compared to the SM. With all these new 
parameters, there are many possible new sources for CP violation. The phenomenology 
of CP violation caused by these new sources is rich and diverse. The effect of these new 
sources of CP violation may surface in the data before any super-particle is discovered. 

Even in the minimal supersymmetric standard model (MSSM), which only augments 
superpartners of known particles in SM, the Higgs sector contains new sources of CP 
violation in its couplings to super-particles. When the fi term in the Higgs superpotential 
and the soft-SUSY-breaking A terms are complex, the tri-boson-couplings between the 
Higgs bosons and the squarks can contain CP violation. In MSSM with the simplest 
universal soft supersymmetry breaking^, there are two new CP violating couplings which 
can be defined to be the phases of fi and A terms in a convention that makes the others 
new couplings real. Therefore, these new sources of CP violation are generic of any 
supersymmetric theories. In addition, they also have been used as one of the leading 
sources of CP violation in a scheme to use MSSM to generate baryon number asymmetry 
of the Universe in electroweak phase transition Therefore, it should be important 
to look for collider phenomenology that can check these mechanism. For example, these 
complex couplings lead to a complex phase in the mixing of stop states. It is the purpose 
of this paper to investigate one consequence of this CP violating source in colliders. 

It is expected that the future colliders are able to produce CP violation signals || [7], 
rOL Of in the sectors of heavy particles. In this article we study the CP asymmetry in the 



Higgs decay into top pairs because the large top or stop coupling to the Higgs particles 
can produce largest effect. 

In MSSM, even with soft breaking terms and R symmetry breaking terms, there is 
no tree level mixing between the scalar and the pseudoscalar bosons. Therefore their 
couplings can be discussed separately. However, the scalar and the pseudoscalar bosons 
mix at one loop, and their effect has to be taken into account as we will show later. 



2 Stop Mixing 

The source of CP violation that we investigate here is due to the mixing in the stop mass 
matrix. We use the convention adopted in Ref. ||. The mass matrix for the stop quarks 
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in the left-right basis is given as 

,2 i ™2 



Mi 



m Q + m t + A t - i m|cos2 / 9 -m t (/xcot/3 + ^) \ 
-m t (/2* cot /? + A t ) + m 2 + A t ~ H m| cos 2/? y ' 



where A^ = | — | sin 0^, and A t ~ fl = — | sin W . The complex phase 5 of the off-diagonal 
elements is the source of CP violation. 



ji* cot (3 + A t = cot (3 + A t \e 



iS 



(2) 



The stop mass eigenstates, ti, t 2 , are related to the left and right stop states by an unitary 
mixing matrix 
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The masses of these eigenstates are given by 

777,2 2 = 1 + m ^ + 2m 2 + ^ 1 _ 4 ^2 ^) m | cos 2/5 =p v^R) , 

i?= ^rrig - + cos 2/3^ + 4m 2 | / u cot (3 + A* t \ 2 
Here we denote t\ as the lighter state. The mixing angle is given as 



tan# 



[m 2 Q — mfj + \m\ cos 2(3 + y/R] 



2m t \fi cot (3 + A* t | 

Strong gluino couplings to stops and tops in the left-right basis is given by 



C § = -\p2g s t R gt R - \p2g s t L gt L . 



(3) 

(4) 
(5) 

(6) 
(7) 



In terms of mass eigenstates, 



cos 9 — P L sin 9e iS )t 1 g - V2g,t(P R sin 9 + P L cos 9e lS )t 2 g, (8) 

where Pl is the left projection \{1 — 7 5 ), and P R is the right projection |(1 + 7 5 ). We 
also need the stop-stop coupling to Z, 



Cz = 
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sin 2 9 w )i{ id" t L - | sin 2 9 w t R id" t 



cos 9w 

After mixing, the Lagrangian for the Z coupling in the t±, t 2 basis is 
9 



R 



(9) 



C 7 = 



cos 9 



w 



(j cos 2 9-1 sin 2 9 W ) t\ id" h + (| sin 2 9-\ sin 2 9 W ) t\ id" t 2 (10) 



+ \ sin(2#) (i\ id" t 2 + t[ id" h 



Note that the last term is real in this phase convention. 
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3 The Higgs couplings to stops 

In MSSM, there is only one pseudoscalar boson., A . The pseudoscalar Higgs coupling to 
the stop squarks is given by the Lagrangian 

' t 



Ca = (t{ 9) T A I .A° (11) 



The matrix T A is given as 



vi _m t ( 2 sin 9 cos 9 Im(i) -z(cos 2 9 A* + sin 2 OA) 
v 2 \ i(cos 2 9A + sin 2 6>A*) -2 sin 9 cos 9 Im(A) 



and A is defined as A = (A t cos /? — //* sin (3)e~ l5 . Note that the nonvanishing of Tfj or T 2 ^ 
is a sure sign of CP violation already (similar to Kl — > 2tt). However, if for some reason 
\i and A t happen to have the same phase, Tj\ and T 2 4 will vanish because in this very 
special case the phase in the stop mass matrix and that in the pseudoscalar couplings can 
be removed simultaneously 

The pseudoscalar Higgs coupling to the top quark is given by the following Lagrangian, 

C\ = JEH- cot pti^tA . (13) 

The neutral scalar Higgs sector is made up two scalar eigenstates, H° and h°. There 
masses are given as 



m2 H,h — \ 



m 2 A + m| ± J (m A + m 2 z ) 2 — Am 2 A m 2 z cos 2(3 



(14) 



Since in MSSM the constraint on the lightest scalar, h, is such that it is too light to decay 
into the top pair, we shall concentrate on the decay of the heavy Higgs boson, H, which 
can decay into the top pair. Our general framework can also be used for the decay of the 
lighter boson, h, of course, if for any reason that it should be heavy enough. The heavy 
Higgs coupling to the stops in the left-right basis is given as 

/ _j^ A coa f a + p\_gn^B (M sin a + a cos a) \ , s 

°~l ^(Asin« + /i*cos«) — cos(a + /3)-^^ J ' 1 ' 

\ 2 sin p v L • ' cos By/ l R v ' / m-\y sm p / 

where the mixing angle a is given in Ref . || . This matrix must then be transformed into 
the stop mass eigenstates. This is accomplished by using the stop mixing matrix, 

T H = U f T^U. (16) 

Its Yukawa coupling is 

„ v qm t sin a _ „ , . 

C Y H = — ttH°. 17 

H 2m w sm(3 K ! 
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4 Helicity calculation of the matrix element 



To get non-zero CP asymmetry in Higgs decays, in addition to CP violating couplings, it 
is necessary to get the absorptive parts from the decay amplitudes in order to overcome 
the constraint from the CPT theorem. We labeled S 1 and P 1 as the absorptive form 
factors of the Higgs or pseudo-Higgs couplings to the top quark. They begin to appear at 
the 1-loop level, unlike their dispersive parts S and P, which can exist at the tree level, 

M = u{p) [{S + iS T )l + i{P + iP 1 )^] v(p') . (18) 

In the Weyl representation, the 7 matrices are given by 





The free spinors of momenta p,p' and helicities A, A' are given by 

where the x' s are two component spinor eigenfunctions a ■ p X\(p) = ^ XX- The u± are 



functions of the energy and momentum of the particles, u± = J E ± \p\. Notice that the 
helicities of tt must match A' = A because of conservation of angular momentum. Our 
normalization of the spinor is u\u\ = v\v\ = IE. The asymmetry between the left and 
right matrix elements is given by 

\M LL \*-\M RR \> 

\Mu\ 2 + \MbrF' 1 J 

The matrix elements are given by 

M LL = ^{-/3t(S + iS I )-t(P + tP 1 )], (20) 
M RR = ^[-PttS + iS^+itP + iP 1 )} , (21) 

with fa — (1 — 4m t /s)2 and s = m? H . The asymmetry can finally be obtained using the 
definition from Eq. fll9|) , 

2(3 t (PS J - P^) 

pl2 + P 2 + ffi S 2 + 02 S I2 ■ ^> 

Since we assume the Higgs boson has definite CP parity at the tree level, the final state 
interactions due to exchanging gluons or gauge bosons in Ref. || are not able to produce 
this CP asymmetry at the one-loop level. However, the rich CP phases in the sector 
of SUSY partners, especially the gluino and the stop, can give rise to large A. For 
scalar boson decay, the second term P 1 S in A gives the leading contribution; while for 
pseudoscalar boson decay, the first term, PS 1 is the leading contribution. 

The polarization asymmetry is Eq.fl2"2"|) can be translated into the lepton energy 
asymmetry || |Tl|, [Tj| in the final semileptonic channel t —>■ hi + v. The energy Eq(£ + ) 
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Figure 1: Triangle diagram via gluino exchange 



distribution of a static t quark decay t — > £ + ub is very simple in the narrow width Tw 
approximation when m& is negligible. 

, j x (l - x )/D if m 2 w /m 2 t <x <1, , . 

/(;Eo) " j otherwise. (23) 

Here we denote the scaling variable Xo = 2E (£ + )/m t and the normalization factor D = 
| — |(miy/m t ) 4 + \{m w /m t ) & . When the t quark is not static, but moves at a speed (3 t 
with helicity L or R, the distribution expression becomes a convolution, 

-'/(1-/30 &Zq ± (s - Sp) 

/./(i+A) 2xl(5i 



fR,L^Pt) = j /fro) 0X2 • (24) 



Here a; = 2E(£ + )/ E t . The kernel above is related to the (lztcos^) polar angular distribu- 
tion. Similar distributions for the i decay is related by CP conjugation at the tree-level. 



Using the polarization asymmetry formula in Eq. (|22|) , we can derive expressions for the 
energy distributions of i~ and £ + : 

N^dN/dx^) = \{l±A)f L {x,(3 t ) + \{lTA)f R {x,p t ) . (25) 

Here distributions are compared at the same energy for the lepton and the anti-lepton at 
the rest frame of the Higgs boson, x(£~) = x{£ + ) = x = &E(l )/Mjj. To prepare a large 
sample for analysis, we only require that each event has at least one prompt anti-lepton 
£ + from the t decay or one prompt lepton £~ from the t decay. 

5 Absorptive parts of 3-point vertices 

We first study the triangle diagram via gluino exchange in Fig. 1.. 
5.1 The Mn stop loop 

When the intermediate state is txti, the Feynman rule gives 

iMu = (-^.) 2 j^ k2 : mika '_^ j! iTnC F , (26) 

where An is defined as 

An = u(p)(P R cos 6 - P L sin 9e iS ) % ^{ + ^ (P L cos 9 - P R sm6e- iS )v(p'). (27) 

q 2 — m| 
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The color factor Cp is |. The absorptive part of the amplitude which is needed for CP 
violation is obtained by cutting across the momentums k and k! . The discontinuity JT3[] 
of the matrix element is 

o 2 T 

Disc(iM u ) = ^-JLfrCp 
8n 



1 — 75 cos 29) + nig sin(2#)(— cos 5 + 27 sin 5) , dQ 



J ' q 2 - nil 47T 

The phase space integration involves the following forms, 

"T^^ = -^(P - pT + ^(p + P'T- (30) 
— m~ 47r 

Multiplying both sides by (p — p')^ the -f/y function can be isolated out because (p + p') ■ 
(p — p') is zero. The Hij function for any intermediate mass and rrij is 

#fT y = 1 + _ # + 4m|/ S )^ . (31) 

The Pij function is given by 



fa = yjl - 2(m 2 + m])/s + (m 2 - m 2 ) 2 /s 2 . (32) 



Notice that when i — j, Pij reduces to = Jl — Amf/s. The function is obtained 
by contracting Eq. (^) with p + p'. 

K ij = -\ J ij( m2 i ~ m ))/ s ■ (33) 

Notice that for matrix elements with both the stop and the anti-stop of the same type, 
the term proportional to (p + p') in Eq. (|30|) does not contribute. 

After this, the imaginary part of the matrix element is needed. The imaginary parts 
are obtained by using the relation 

Disc(jM) = 2iu(p) \S T 1 + P 7 Z7 5 1 v(p') . (34) 



s i = 9]TuPi Cf ( +m 7 sin ( 2 g) C os5 + 2m t H u ) , 

167TS 
167TS 

5.2 The M22 stop loop 

Similarly, we obtain results for the intermediate state £2^2- The matrix element is given 
by 
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where A22 is given by 

M 22 = u(p) (P R sin 9 + P L cos fle ig ) ' ( ( + m £ (P L sin 9 + P R cos 6e~ t5 )v(p') . (37) 

g z — mj 

After integrating the phase space of the intermediate state in the cut diagram, the form 
factors are 

g 2 T 22 p 2 , 



22 



IQlTS 



p i 9 2 S T 22 (3 2 

r 22 — 



(—m~ g J 22 sin(2#) cos 5 + 2m t H 22 ) 
Cf (+mgJ 22 sin(26 l ) sin5) . 



(38) 



167TS 

5.3 The A^i2 stop loop 

The amplitude involving the intermediate state t±t 2 is given by 

d 4 q 



For t 2 tx, it is 
iM 21 

M 21 



2fglT 12 C F [ M l2 1 1 

J k z — mf k' 2 — m 2 {2n 



u(p)(P R cos9 -P L sm9e iS ) 



J+ m g) 
q 2 — m~ 



(P L sm9 + P R cos9e- iS )v(p'). (39) 



(27T) 4 ' 



(40) 



u(p)(P R sin + P L cos 9e i5 ) ( 4^-^4(P L cos - P fi sin 0e~ i5 W) . 

q 2 — m~ 



After integrating over the intermediate phase space, we add up the absorptive parts to 
give 



°12+21 



P J 

-'12+21 



8tts 
8?rs 



C F Re m~ g T 2l (cos 2 #e i<5 - sin 2 0e _ld ) J 12 , 



C F Im 



m- g T 2 i (sin 



cos 2 fle* 5 ) 



+m t sin(26 , )T 2 i(m2 — m\)/s 



J 12 ■ 



(41) 



(42) 



6 Absorptive parts of 2-point vertices 

We study the bubble loops which involve only the it pair. 



6.1 Z diagrams 

Z diagrams that contain t\ti or t 2 t 2 are identical zero because of the phase space inte- 
gration^ This point will become obvious from the result of the mixed intermediate states 
t\i 2 or t 2 t\. The A4i2 matrix element is given below, 



iMv> = 



g 2 T 12 N c L 
4 cos 2 9w 



sin(2£) J Uz-^ 



1 



d A k 



m 2 z k 2 -m\ k' 2 -m 2 2 (2tt) 4 ' 
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Figure 2: Z exchange diagram 



A A H° 




Figure 3: Higgs mixing diagram 

where Mz is given as 

M z = u(p)Y (| - | sin 2 6 W - i 7 5 ) v(p') (g, v - y,/m|) (k - k') u . (43) 

The M.21 matrix element is very similar to the above matrix element, with the substitution 
ofT 12 byT* 2 , 



p j 9 2 N C m t 12 sin(2fl) m\ - mf 

Pl ™- 647T mlcos^^ ^ Im ™, (11) 



where the color factor Nq = 3. This graph will contribute only to CP violation of the 
scalar Higgs decay. One may think that without gluino couplings in the graph, one 
should be able to rotate away the CP violating phase in scalar coupling Tf^. However, 
such rotation will produce a complex phase in t[t 2 Z coupling in Eq.(lO). 

6.2 A°-H° Mixing 

The stop bubble loop induces A°-H° Mixing. We study its absorptive part which con- 
tributes to the CP violation. In the heavy Higgs mass limit of MSSM, m^o and m#o are 
quite close to each other based on the tree- level mass relation in Eq.(|lJ]). However, it is 
known that there is large higher order correction to the tree-level mass relation. Thus 
in our following study we allow the masses m A o and m H o to vary independently, not re- 
stricted by the tree-level formula. The matrix element for the pseudoscalar Higgs turning 
into a stop pair, then becoming the heavy Higgs, and finally decaying into a top pair, is 
given as 

igrrn sin a y f T,f T£ *A/(2tt) 4 

%JVi 2m w sm(3 U[P)V[P) j-J m\ - m\ (I + q)* - m 2 (I - qf - m) ' ^ 

Making the same cut as that in the Z loop digram, we obtain the imaginary part of the 
form factor, 

S'(A° it) = _^!^g^ EA T / T 4° . (46) 
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A similar expression is derived for the heavy scalar Higgs decay, 

P\H« - A - it) = --f^- cot (3 £ fa, J* \ . (47) 

32^^ ^ m 2 H0 - m 2 A0 

7 Physical and Numerical Analyses 

Before we plunge into the numerical analysis, it is interesting to check the limit in which 
the two stop states are accidentally degenerate. In that case, (A4?)x2 = («M|)2i = 0, 
(•Mf)n = CMf)22- Therefore fi* cot (3 + A t = 0, and fi* and A t should have the same 
phase which can still serve as the source of CP violation. In that case, 9 and 5 in U in 
Eqs.(3),(12),(16) and in the definition of A, should be set to zero. 

Thus in this limit, the stop loops do not contribute to M.\\ and M.22 in the pseudoscalar 
case, because Tj\ = T 2 ^ = 0. However they still give rise to CP violation in .M12, M.21 
because = (— im t /v2)(A^ cos/? — /isin/3). One may attempt to absorb this phase by 
rotating the phase of, say, the right stop, however such rotation will lead to complex gluino- 
top-stop couplings which cannot be rotated away because of the nonvanishing gluino mass. 
From this, it is easy to understand why a factor of gluino mass has to appear in Eq.(41) 
for S[ 2+2 i- Similarly, for the scalar Higgs decay in the degenerate stop limit, the stop 
loops still produce no CP violating effect in A4u and M.22, because sin6> = and only 
the term proportional to the gluino mass in P/2+21 contributes as reflected in Eq.(42). 

It is also straightforward to note that in the degenerate limit, the contributions of both 
H°-Z° and A°-H° bubble graphs vanish. In the H°-Z° case, the phase of the scalar Higgs, 
H, coupling as well as that of the stop mixing can be rotated away simultaneously (into 
the gluino couplings) without affecting the Z coupling and this is reflected in m\ — m\ = 
factor in Eq.(44). For A°—H°, the phase of pseudoscalar coupling as well as that of the 
scalar coupling can be rotated away simultaneously also and this is reflected in 

EfriTfTf = P 12 (Ti° Tl f + Ti°T&°) = , 

ij 

in this particular limit. 

To illustrate our result numerically, in the following, we set the parameters so that 
only the lighter stops states t\ti are light enough to be on-shell for simplicity. In such 
a scenario, only some of the above contributions are available. In Fig. 4, we show the 
mass rrii of the lightest stop versus tan f3. The best current limit of the lowest bound on 



lightest stop mass from LEP is about 95 GeV [|14| , and this means that tan (3 < 3 is not 
allowed for the case mg = vrtu = 300GeV, [i = 500 GeV, A t = 500^ GeV. If one wishes 
to study the possibility of a much heavier Higgs which can decay to all channels of stops, 
the remaining diagrams can be easily incorporated into the numerical analysis. 

7.1 Pseudoscalar-Higgs Decay 

In the model of our study, the A remains its status as a pseudoscalar boson at the tree 
level. CP violation in the pseudoscalar Higgs decay into top pairs occurs starting at the 
one-loop level. The leading contribution requires induced scalar form factor S 1 which, 
as we have shown, can be obtained from the absorptive part due to the intermediate tt 
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state. Notice that there is no Z loop contribution to S 1 in the Higgs decay. Fig. 5 show 
the asymmetry for the pseudoscalar Higgs decay as defined by (p2]). Fig. 6 shows the 
branching ratios of the pseudoscalar-Higgs decay to top pairs, bottom pairs, and stop 
pairs. For small tan (3 the decay channel is mostly top pairs. 

7.2 Higgs Decay 

In the Higgs decay, the CP violation is caused by terms proportional to the P 1 form 
factors. The Z diagrams can contribute in principle if not disallowed by the kinematics. 
As stated before it does not contribute in our illustration because we assume a heavy fa- 
Fig. 7 shows the CP asymmetry of the Higgs decay. The branching ratios for the Higgs 
to decay into tops, bottoms, stops, W's, and Z's are given in Fig. 8. 



8 Conclusion 

The complex mixing among the stop sector can produce CP asymmetry at the level of a 
few percent in the final products of polarized tt states from the Higgs boson decay. Such 
asymmetry can be measured in the energy spectra of the final leptons. Unlike the usual 
two-Higgs-doublet model, the CP violation does not require the mixing among A and 
H° states at the tree level. 

This work was supported in parts by National Science Council of R.O.C., and by U.S. 
Department of Energy (Grant No. DE-FG02-84ER40173). 
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Fig. 4 Light stop mass m x in GeV versus tan j3 for the case rriQ = ra v = 250, 300, 350 GeV, 
/j = 500 GeV, A t = 500ex GeV. 
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5 Asymmetry of pseudoscalar Higgs decay itlq = 300 Gev, mjj = 300 GeV, \x = 
500 GeV, A t = 500ex GeV, m A = 400 GeV, m H = 420 GeV. 
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6 Branching ratios for pseudoscalar Higgs decay, m Q = 300 Gev, = 300 GeV, 
H = 500 GeV, A t = 500ex GeV, m A = 400 GeV, m H = 420 GeV. 
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7 Asymmetry of heavy Higgs decay, vtlq = 300 Gev, rriu = 300 GeV, fj, = 500 GeV, 
A t = 500ex GeV, m A = 400 GeV,m H = 420 GeV. 
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8 Branching ratios for heavy Higgs decay, uiq = 300 Gev, mu = 300 GeV, \i 
500 GeV, A t = 500ex GeV, m A = 400 GeV, m H = 420 GeV. 
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